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1. Introduction 

The study of neutron star properties such as their masses and radii has been a topic of high 
interest over the last decades. Very recently, the measurement of the Shapiro delay for the mil- 
lisecond pulsar PSR J1614-2230 has provided a new reliable limit for the highest known pulsar 
mass of (1.97 ±0.04) M e [p. The value for the maximum mass is of major interest not only for 
astrophysics but also for nuclear physics as it is closely connected to the properties of dense nuclear 
matter, especially the stiffness of the nuclear matter equation of state (EoS). Nuclear matter at finite 
baryon number densities is also subject of studies at heavy ion facilities such as Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven, the future Facility for Antiproton and Ion Research (FAIR) 
at GSI, Darmstadt, and the Nuclotron-based Ion Collider facility (NICA) in Dubna, as well as the 
Facility for Rare Isotope Beams at the Michigan State University in East Lansing, USA. 

This work aims at studying the impact on the properties of compact stars of a soft nuclear 
matter EoS, as obtained from the KaoS experiment at GSI, Darmstadt, for baryon densities up 
to two-three times nuclear matter saturation density uq ~ 0.17 fm~ 3 . The K + meson production 
in nuclear collisions at subthreshold energies is used to study the level of compression of dense 
matter, which, in turn, is controlled by the stiffness of nuclear matter through the nucleon potential 
Ujv- The more attractive U^ is, the higher is the produced K + meson abundance, which therefore 
can serve as a probe for the stiffness of nuclear matter [Q] . 

This paper is organized as follows. First, we review the results on K + meson production from 
the KaoS experiment and the need of a soft EoS to explain the observed data. Next, we study 
the relevance of the KaoS results for low mass neutron stars and their connection to the nuclear 
symmetry energy. This question arises since light neutron stars can have interior densities in the 
same range as the ones probed by KaoS. The last part of the paper is focused on the determination of 
the upper limit on neutron star maximum masses by using the information coming from the KaoS 
experiment and causality constraints. We conclude the paper with a discussion of the obtained 
limits for the highest possible neutron star mass with respect to neutron star observations, such as 
the PSR J 1748 -2021 pulsar 

2. A soft nuclear equation of state from KaoS experiment 

The measurements on K + meson production in nuclear collisions at subthreshold energies 
were performed with the Kaon Spectrometer (KaoS) at GSI. The beam energy dependence of the 
K + multiplicity ratio from Au+Au and C+C collisions at subthreshold energies of 0.8 to 1.5 GeV 
per nucleon was introduced as a sensitive and robust probe for the stiffness of nuclear matter [Q]. 

In order to describe the experimental results, IQMD and RQMD (Isospin and Relativistic 
Quantum Molecular Dynamics) transport model calculations were done [||, g], applying a Skyrme 
type Un with two- and three-body forces. The detailed analysis showed that the K + multiplicity ra- 
tio from Au+Au and C+C collisions is rather insensitive to cross-sections, momentum dependence 
of in-medium potentials and transport models used while the experimental uncertainties cancel out 
in the double ratio [Q]. The transport calculations consistently demonstrate that the beam energy 
dependence of the kaon multiplicity ratio is described best by the nucleon potential with small 
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repulsion. When applied to infinite isospin symmetric nuclear matter this conesponds to a soft 
nuclear EoS characterized by a stiffness parameter K = 9^ |„ = 200 Me V [||]. 

The measured attractive nucleon potential and the corresponding soft nuclear EoS [Q, ||, ^] 
have to be tested on their compatibility with neutron star observations. In this work we want to 
analyze the implications of the results from the KaoS collaboration on compact stars. 



3. Radii and moments of inertia of low mass neutron stars: 
density dependence of the nuclear symmetry energy 

The neutron star radii are a focus of major attention [|7| ||, ^, [l(], 11, 12, 13, [14|, |15|] due to 
their implications for nuclear matter properties. Also moments of inertia seem to have captured the 
interest of the astrophysics community. Guillemot et al. have presented a new technique to deduce 
limits on the moment of inertia of gamma-ray pulsars This analysis is based on the efficiency 
of the gamma ray emission from gamma-ray pulsars with respect to the total energy loss obtained 
from the pulsar spin-down. Due to the fact that the energy emitted in gamma-rays cannot be larger 
than the total emitted energy, one can obtain a lower limit on the pulsar's moment of inertia. Also 



moments of inertia have been addressed in connection with double pulsars p7|]. Such studies could 
be substantially refined in the future with the Large Observatory For X-ray Timing (LOFT) [18]. 

Neutron star radii and moments of inertia are strongly connected to the symmetry energy [|T^] . 
For our study, light neutron stars are of special interest, because their central densities can be in 
the same range as the ones probed by the KaoS experiment. Since the latter provides information 
on the stiffness of nuclear matter, the symmetry energy is the remaining uncertainty of the nuclear 
EoSs. 

For a consistency, we choose a Skyrme type EoS with a nucleon potential similar to the one 
which was used in the analysis of the KaoS data. The energy per baryon is written as [ |20| ] : 

| = m n (1 -Y p ) +m p Y p +E J +D " 

+ (l-2F / ,) 2 [(2f-l)£ («5-f(«))+5o« r ], (3.1) 

whereas u = n/no is the baryon number density and Eq is the average kinetic energy of nuclear 
matter at no and a proton fraction of Y p = 0.5. Two- and three-body forces, described by the terms 
B and D, together with a, are fitted to reproduce the binding energy per baryon E(no,Y p = 0.5) = 
— 16MeV, the stiffness parameter K, and the saturation density uq. While So is varied between 
28 MeV and 32 MeV, its density dependence is chosen as a power law with u r where 7 = 0.5 — 1.1 

1,0 in, p. 

Eq. ( |3~l| ) is used to describe matter in the neutron star core, while we incorporate an inner and 
outer crust j f24] , [25| ]. We first study non-rotating neutron stars with M = 1.25 M , in accord with 
the lightest pulsar masses [|27|]. The radii and moments of inertia in dependence of K are shown 
in Fig. [j] for different symmetry energy configurations. It can be seen that both largely depend 
on the density dependence of the symmetry energy given by y, i.e. L = 3no(dS(rib) /dii},)\no. At 
K ~ 200 MeV, stiff and soft symmetry energy configurations lead to a difference in the neutron 
star radius and moment of inertia of AT? ~ (1 — 1 .5) km and AI ~ (2 • 10 43 — 2.5 • 10 43 ) g cm 2 , 
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Figure 1: Radii R and moments of inertia / of non-rotating neutron star configurations with 1.25 M Q as 
function of different slopes of the symmetry energy L. Central densities are given in units of the nuclear 
matter saturation density Hq. Taken from $2M. 
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Figure 2: Radii and moments of inertia of light neutron stars for Kq ~ 200MeV as function of their mass 
for different symmetry energy setups as in previous figure. Central densities are given in units of the nuclear 
matter saturation density «o- Taken from [E8fl. 



respectively. The central densities of the corresponding stars are in the range of < 3.4 no, similar 
to the density region explored by the KaoS collaboration. 

In Fig. ^, we study the radii and moments of inertia of light neutron stars with masses in the 
range of (1.1 — 1.6) M for Kq ~ 200 MeV and with equal configuration for the symmetry energy 
as in Fig. [l| It can be seen that especially neutron stars with M < 1.3 M have central densities 
which are low enough so that the entire neutron star interior can be described by an EoS which is 
probed by the KaoS experiment. With the isospin symmetric part of the nuclear EoS determined 
by the latter, radius and moment of inertia measurements from light neutron stars could thereby 
distinguish between a soft and a stiff behavior of the symmetry energy. 

4. Maximum mass for neutron stars 

The most massive stable neutron star configuration depends on the stiffness of nuclear matter. 
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Figure 3: The nucleon potential as a function of the baryon number density for different models. The 
vertical thin lines indicate different densities for the onset of the stiffest EoS. Taken from 



The softer the nuclear EoS, the lower is the maximum mass which can be reached before the star 
becomes unstable against collapse to a black hole. 



Ref . [ 26 ] introduced the idea to use known properties of hadronic matter up to a critical density 



n cr i t . At higher densities, the low density EoS is smoothly connected to the stiffest possible EoS 
allowed by causality. Previous studies showed that the highest possible mass can be expressed as 



M high = 4AM Q (e crit /£o) 



-1/2 



(4.1) 



with Bait being the energy density taken at the transition density n cr [ t from the known low density 
EoS to the stiffest one. In Refs. [29, 31, 32] an extrapolation of measurements on nuclei at n ~ n 
to higher densities were made and resulted in a maximum mass of M/„- ? /, ~ 2.9 M for e cr j t = 2 £ - 
KaoS data directly tests the EoS for densities n ~ (2 — 3) no and justifies the choice of n cr i t in the 
same range. Consequently, we will calculate the highest possible masses with the restriction on 
n cr it given by the KaoS data. 
For 

ftcrit up to (2 — 3) wo the nucleon potential should fulfill the results from the analysis of the 
KaoS data. For Skyrme type models this restriction corresponds to K < 200 Me V. As seen in Fig. ||, 



we use two representative approaches for Skyrme models, BSk8 and SLy4 [33, 34]. For Relativistic 
Mean Field (RMF) EoSs, the stiffness of high density nuclear matter is determined by the nucleon 



effective mass m* at uq [35]. Consequently, the corresponding Um is chosen by varying m* for 



a given compression modulus Kq, so as to obtain a nucleon potential which is similar to or even 
more attractive than the Skyrme parametrization within the density limits. A more attractive Un at 
supra-saturation densities allows a higher compression of matter for the same bombarding energy, 
which enhances multiple-step processes in subthreshold kaon production and is in line with KaoS 
results. Fig. [?]also shows that the TM1 [36] parametrization as well as the Brueckner-Hartree-Fock 
approximation (BHF) [J37], H H, EH fulfm this requirement. Schemes, such as BSk8 [||], SLy4 
[ 34 ] or NL3 pi] ] , produce nucleon potentials which are more repulsive than the Skyrme benchmark 
in the density region of interest. 

Fig. [I| shows the highest possible neutron star masses at n cr j t ~ (2 — 3) «o for the discussed 
EoSs. Note that the maximum mass does not depend on the asymmetry at low densities but on the 
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Figure 4: Highest possible masses for compact stars for different EoS as a function of the critical density 
for the onset of the stiffest possible EoS. Taken from |[28||. 



critical density. The higher n cr j t is, the later is the onset of the stiffest EoS in the star's interior. 
Consequently, less mass is supported and the value for the highest mass decreases. For n cr n ~ 2?io, 
the star is dominated by the stiffest EoS and therefore reaches masses of up to 3M . Smaller 
maximum masses of Mugh ~ 2.4 M Q are obtained for the upper limit of n cr j t ~ 3 no. However, 
this lower bound depends on the EoS studied. Moreover, it can be observed that a lower bound of 
2.8 M as given by the NL3 EoS is not compatible with constraints from KaoS data. From Fig. |] 
we can conclude that a pulsar mass of 2.7 M Q as found for PSR J1748-2021 by [||] is marginally 
compatible with a soft EoS and requires a prompt transition from a soft EoS to the stiffest possible 
for a density around (2.2 — 2.5) no- 

5. Summary 

We have studied the properties of low mass neutron stars, such as radius and moments of 
inertia, using K + multiplicities from heavy-ion collisions at GSI. Light mass neutron stars can 
have interior densities in the same range as the ones probed by KaoS experiment. Thus, these 
experimental results serve to constrain the stiffness of the nuclear EoS for low mass stars. The radii 
and moments of inertia of low mass neutron stars are then sensitive to the density dependence of 
the symmetry energy, which is the remaining uncertainty of the nuclear EoSs. We find that light 
neutron stars withM < 1.25 M are well suited for future radius measurements, such as the LOFT 
mission. 

Moreover, we apply the KaoS results up to densities of 2 no < n cr j t < 3 «o and introduce the 
stiffest possible EoS for n > « cnf to calculate the highest allowed maximum mass. In this way, 
we test whether a soft nuclear EoS up to ~ 3 no is compatible with tentatively massive neutron 
stars such as PSR J1748-2021. The KaoS results confirm the previous theoretical estimation for 
the highest possible neutron star mass of ~ 3M , estimation based on an extrapolation to supra- 
saturation densities of results for nuclei at n ~ no. A pulsar mass of 2.7 M is not excluded by the 
KaoS data, but requires the onset of the stiffest possible EoS already at n cr j t ~ 2.2 — 2.5 no- Future 
experiments, such as FAIR at GSI, aim at probing matter at densities beyond 3 no and, therefore, 
giving further constraints on maximum masses. 
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